Unraveling the "Pressure-Effect" in Nucleation 
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The influence of the pressure of a chemically inert carrier-gas on the nucleation rate is one of the biggest 
puzzles in the research of gas-liquid nucleation. Different experiments can show a positive effect, a negative 
effect, or no effect at all. The same experiment may show both trends for the same substance depending on 
temperature, or for different substances at the same temperature. We show how this ambiguous effect naturally 
arises from the competition of two contributions: nonisothermal effects and pressure-volume work. Our model 
clarifies seemingly contradictory experimental results and quantifies the variation of the nucleation ability of a 
substance in the presence of an ambient gas. Our findings are corroborated by results from molecular dynamics 
simulation and might have important implications since nucleation in experiments, technical applications, and 
nature practically always occurs in the presence of an ambient gas. 

PACS numbers: 82.60.Nh, 64.6().qe, 64.60.Q- 



A gas phase can be supersaturated considerably beyond 
its equilibrium condensation point before liquid drops form 
spontaneously. The pathway of the phase transition is blocked 
because microscopic droplets are thermodynamically less fa- 
vorable than the bulk vapor. Therefore, the transition can 
only be initiated by rare fluctuations exceeding a critical size, 
called the critical nucleus. The formation of such a critical nu- 
cleus is the limiting step in the transition and its frequency of 
occurence is called the nucleation rate. Nucleation is behind 
most phase transition and plays a crucial role in atmospheric 
processes such as the formation of aerosols or the condensa- 
tion of water vapor into clouds [1, 2]. An accurate experi- 
mental evaluation of the nucleation of atmospherically rele- 
vant substances and its correct theoretical prediction are es- 
sential for a better understanding of climate change and are 
the subject of intense investigations [3]. 

Nucleation is highly sensitive to small changes in the state 
variables describing the system, most notably to temperature. 
However, condensation is inevitably connected with the re- 
lease of latent heat. In experiments this latent heat is removed 
by the presence of a large background of caiTier-gas. This 
carrier-gas should be noncondensing and chemically inert and 
should have no influence on the nucleation except serving as 
the desired heat bath. Nevertheless, many experimental re- 
sults suggest that there can be an influence on the nucleation 
rate that can span some orders of magnitude, depending on 
the pressure and type of the carrier-gas (see Refs. (3, '5] and 
references therein). Unfortunately, the overall picture is far 
from clear Comparisons of experimental findings cover all 
possibilities: no effect, increase, or decrease of the nucleation 
rate with carrier-gas pressure. There were some concerns if 
the pressure-effect is not just an experimental artifact, which 
cannot be completely ruled out in all cases. However, continu- 
ous improvements in the experimental setups rather confirmed 
than remedied this undesireable and elusive effect. 

Following the experimental evidence, theoretical works 
have tried to explain the pressure-effect with varying suc- 



cess. These approaches include changes to classical nucle- 
ation theory (CNT) accounting for nonideal behavior of the 
vapor and/or caiTier-gas |6], treating the carrier-gas as not 
fully inert or the problem as binary nucleation |7, 8, 9] or an- 
alyzing its influence on cluster stability and the impingment 
rate 1 10], and many more. There remain many contradictions 
and differences both in the direction and in the magnitude of 
the effect. These are complicated by the discrepancies be- 
tween experimental results themselves and in the comparison 
of experiment and theory. 

Here we present a simple yet physically very appealing 
model that resolves many of the apparent contradictions of 
the pressure-effect. We take a deliberate step back and incor- 
porate the presence of a carrier-gas into classical theory in a 
most natural manner that accounts for the two primal contri- 
butions of the carrier-gas: the efficiency of thermalization and 
the additional work that a cluster has to spend for growing in 
its presence. These contributions have opposite trends and we 
show how this may be responsible for the existence of appar- 
ently contradictory results. 

Nucleation theory [11] usually aims at calculating isother- 
mal nucleation rates, taking the constant value of the tempera- 
ture and the idea of a non-influential carrier-gas for granted. In 
CNT, the work of formation of a droplet of size n at constant 
pressure p and temperature T is L 1 l.l 



AG(n) = -nh.pi + s\Yn 
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(1) 



i.e., the combination a volume term related to the difference 
in chemical potential A/i between the vapor and the liquid and 
a surface term needed to build the liquid interface A — s\n^l^ 

with the surface tension y {s\ — {36kv^) is the surface area 
per monomer and v/ is the volume per molecule in the bulk 
liquid). The vapor pressure p is assumed as ideal and the liq- 
uid cluster is considered as an incompressible spherical drop 
with a sharp interface and bulk liquid properties. The free en- 
ergy has a maximum at the critical size n* and its height AG* 
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AG* = 



16;r v?7^ 
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The isothermal steady-state nucleation rate then is 
/cNT=^exp(AG*AB7'), 



(2) 



(3) 



where K is a kinetic prefactor In CNT, the difference in the 
chemical potential is given by 



Aji =kBTlnS-vi (p-peq) 



(4) 



Here, S — p/ p^q is the supersaturation and p^q is the equi- 
librium vapor pressure. The second term in Eq. is aris- 
ing from the pressure-volume work the liquid drop has to 
perform against the ambient vapor pressure ifisll . It is typ- 
ically small and hence commonly neglected but we keep it 
here for clarity. We now naturally account for the presence 
of an ideal carrier-gas by noting that the cluster must also 
perform pV-woik against the ambient carrier-gas pressure, 
Wc — pcV{n) = m'ipc LMJ. The work of formation including 
pV-work now reads 



AGpv (n) = -nA^eff + ii7n 



2/3 



(5) 



where we have cast the pV-contributions into an "effective 
chemical potential" 



AfieS^^kBTlnS-Vi (p+Pc-Peq) 



(6) 



Replacing A/x for A/Xeff in Eq. (|2|l, the barrier becomes 
AGpY = (167r/3) {yjy^ /AjXeff^) and the pV-corrected rate fol- 



lows directly from Eq. ([3]) as 



Jpv = K exp {AG*y/kzT 



(7) 



The main designated role of a carrier-gas is to keep T con- 
stant. However, this thermalization might not be perfect and 
nucleation then happens under nonisothermal conditions. The 
classical work of Feder et al. offers analytical expressions to 
quantify this influence on the nucleation rate lll5lll6il . Phys- 
ically, it is controlled by the competition between the energy 
increase due to latent heat and the energy removal through 
elastic collisions with vapor and carrier-gas molecules. The 
parameter 
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ksT dA{n) 



dn 



(8) 



quantifies the energy released when a monomer is added to a 
cluster, which is the latent heat h per molecule (corrected by 
a small amount k^T /2) minus the energy spent on increas- 
ing the surface area A(n) against the surface tension. The 
mean squared energy fluctuation removed by collisions with 
impinging molecules is 



(9) 
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Figure 1 : Deviation of the nucleation rate from CNT due to the full 
pressure- effect, Eq. as a function of the ratio of carrier-gas over 
vapor molecules (solid) arising from the two contributions of non- 
isothermal effects (dash, Eq. dlOt ) and pV-work (dash-dot, Eq. 0) 
for Argon at 50 K and 5 = 869. 



for an ideal monatomic vapor and carrier-gas fT^ (m is the 
molecular mass, the number of molecules of the condens- 
able). Eq. ^ indicates that a large number of carrier-gas 
molecules A^c and a light (small molecular mass mc) carrier- 
gas are most effective for a good and fast thermalization. Fi- 
nally, the influence of nonisothermal effects on the steady- 
state nucleation rate is determined by the combination of q 
and b: 



Aionisoth. 



isoth. 



(10) 



We now obtain the full pressure-effect (PE) of the carrier- 
gas on the nucleation rate by combining Eqs. (fTOl i and tak- 
ing the latter as the isothermal but pV-corrected rate: 



JpE — 



b^ 



-1 



(11) 



Obviously, the magnitude of the rate calculated by Eq. (fTTI) 
depends on the estimate of the underlying CNT, which of- 
ten can be off by many orders of magnitude {I?, 18 1. In our 
case however, we are only interested in the deviations arising 
from the pressure-effect. By renormalizing Eq. (fTTI) we get 
a reasonable^ and mostly model-independent estimate of this 
deviation il2ll : 



J, 



pV 



JpE _ I 

Jem b^ + q^ 7cNT ' 



(12) 



Fig. [T] shows the change of the rate as a function of the ra- 
tio of carrier-gas to vapor molecules, Nc/N, which for perfect 
gases is the same as the ratio of carrier-gas over vapor pres- 
sure, pc/ p- The calculation was performed for argon at 50 K 
and a supersaturation of 5 = 869 |19] and the latent heat h 
was calculated via the Clausius-Clapeyron relation. Fig.[T]also 
shows the individual contributions coming from nonisother- 
mal corrections, which always have a positive effect on the 
rate, and pV-work, whose effect is always negative. The com- 
petition of both terms first leads to an increase in the rate with 
increasing carrier-gas pressure due to better thermalization. 
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Figure 2: Comparison of MD results with theoretical prediction of 
the pressure-effect, Eq. l ll2b . 



Then, the penalty of pV-work takes over and leads to an over- 
all decrease of the rate. 

We have performed standard molecular dynamics (MD) 
simulations of Lennard- Jones (LJ) argon nucleation at 50 K 
(Nav = 343, S = 869) and 80.7 K (A^Ar = 512, 5 10.5) with 
LJ-helium as carrier-gas in order to verify the predictions of 
our model. Details of the simulation are similar to Ref. (idn . 
The rates were analyzed using a method based on mean first- 
passage times 1,20,1 . The studied A^c/A^-ratios range from 1 
to 20. We have renormalized the simulation data to the first 
point of each series because we are only interested in the de- 
viation arising from an increase of the amount of carrier-gas. 
Fig. (in shows the simulation result together with the predic- 
tion of Eq. ( fT2] i. which is also renormalized by the same fixed 
value. The agreement is quite remarkable. At 50 K, we only 
observed a slight increase in the rate with increasing amount 
of carrier-gas. On the other hand, the rate already starts to 
drop at a much smaller ratio of Nc/N = 4 at 80.7 K, again 
following the theoretical prediction remarkably close. 

We now take a closer look at the influence of temperature 
on the pressure-effect. Fig. ^ again shows the pressure-effect 
for argon, this time for three different nucleation tempera- 
tures. We took the supersaturations at each T to correspond 
to approx. the same rate from an earlier work |il9J. The non- 
isothermal effects, Eq. ( fTOl i. are only weakly depending on the 
supersaturation of the system. The pV-work term, however, 
gets more and more pronounced the higher the temperature. 
The reason for this is the huge change in the equilibrium vapor 
pressure by almost three orders of magnitude with increasing 
T. Hence, the same pc/p-ratio corresponds to a much higher 
total pressure in the system and the pV-work contribution is 
taking over earlier than at lower T. We have also framed a 
region of ratios that can be encountered experimentally UtIi . 
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Figure 3: pressure-effect for argon at three different temperatures, 
50 K (dash-dot), 65 K (solid), and 80.724 K (dash) as a function of 
the ratio of carrier-gas to vapor. The supersaturations at each tem- 
perature correspond to a base rate of about 10^^ cm^-'s^' 1 19]. The 
vertical lines frame a region of A'c/A'-ratios typically encountered in 
experiments. 



In a typical experiment flS to study the pressure-effect, one 
would fix T and S and vary the total pressure by increasing the 
amount of carrier-gas. We can do the same in Fig. (O, going 
from 3:1 to 20:1. At 50 K we observe a mild increase of the 
rate of about a factor of 2. There is practically no change in the 
rate at 65 K, at least none that would be detectable by available 
experimental techniques. Finally, at 80.7 K, we only observe 
a strong decrease of the rate of up to 2 orders of magnitude. 
Thus, we understand clearly how it is possible to observe only 
a positive, a negative, or no effect at all for the same substance. 
Both the sign at any given pc/ p-ratio and the magnitude of the 
effect will strongly depend on the equilibrium vapor pressure 
(thus on the substance and temperature) and the experimental 
window of accessible pc/ p-ratios. 

Finally, we discuss some of the generic conclusions that 
can be drawn from our model with respect to different ex- 
perimental situations. Since experiments greatly vary in the 
preparation, parameters, rate window, evaluation method etc., 
a careful evaluation requires a comprehensive analysis that we 
leave for a future work. In general, we can distinguish be- 
tween two different limiting behaviors of the pressure-effect, 
depending on the relative influence of the pV-woik term and, 
therefore, the total pressure ptot ~ P + Pc- If ^^f- *C In 5, the 
contribution of pV-work can be neglected and the only no- 
ticeable contribution comes from nonisothermal effects. In 
that case, the possible influence of the pressure-effect follows 
the dashed line in Fig.[T]and the only effect we could observe 
in an experiment would be a slight increase in the rate with 
increasing amount of carrier-gas. Note also that we would not 
observe any pressure-effect at all if we already started from 
high molar fractions of carrier-gas. Interestingly, these are the 
conditions that we typically find in nucleation pulse chambers 
112 ill , where the total pressure varies only very little around 
1 bar and where carrier-gas fractions are high. Our analysis 
thus justifies why no effect is found in these experiments. In 
the opposite limit, when pc/p is large and -^y- ^ InS, the 



system is perfectly thermalized but the dominating pV-woik 
leads to a significant overall decrease of the rate. Moreover, 
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we can make an estimate of the changes in the critical size in 
terms of the CNT estimate 111 ill . 



32k vjf 



(13) 



again replacing the chemical potential by the effective one, 
Eq. (|6]l. Eq. (fTsT i shows that a further increase of the carrier- 
gas pressure (or similarily ptot) would also lead to an increase 
in n*. In an experiment, this might be noticeable by an in- 
crease in the slope of a nucleation rate isotherm, which is 
connected to the critical size via the nucleation theorem iflill . 
However, care must be taken in that kind of analysis because 
the nucleation theorem in principle requires that the carrier- 
gas pressure pc is also fixed. 

In summary, we have presented a simple model that is able 
to shed light on one of the biggest puzzles in current research 
of vapor-liquid nucleation: the "pressure-effect". The inves- 
tigation of this effect is often entangled between experimen- 
tal uncertainties, wildly different experimental conditions and 
procedures, and theoretical interventions on different stages of 
the modeling of nucleation and growth. We have presented a 
physically very appealing way to disentagle most of these am- 
biguities by properly incorporating the presence of a carrier- 
gas into CNT. Simulation results corroborate the validity of 
the model quite impressively. Still, we cannot discard the in- 
fluence of other factors on the observed pressure-effect. Nev- 
ertheless, these factors, if applicable, can be added easily to 
our model as secondary contributions to the more fundamen- 
tal and inevitable physical roles that the presence of an am- 
bient carrier-gas plays in nucleation and which are accounted 
for in our model. For example, it is straightforward to include 
nonidealities of the vapor and carrier-gas as well as the com- 
pressibility of the liquid, which surely will play a more promi- 
nent role the higher the pressure. Another possible influence 
we have neglected here is a change in the kinetic prefactor. 
Finally, we deliberately separated the influence of a truly inert 
carrier-gas pressure from other effects such as binary nucle- 
ation or surface adsorption. It is somehow misleading to in- 
clude these under the same "pressure-effect" tag because even 
though the strength of them may depend on pressure, their ori- 
gin certainly is not the pressure of the carrier-gas. In any case, 
the insights provided by the model presented in this letter will 



undoubtedly be very helpful to quantify and remove the influ- 
ence of the carrier-gas on experiments. This opens the door to 
a more accurate evaluation of nucleation rates, which has im- 
portant implications on many atmospheric and technological 
processes. 
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